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Abstract 

The CP conserving amplitudes for the decays K — > 3ir are calculated in Chiral 
Perturbation Theory at the next-to- leading order. We present the expressions in a 
compact form with single parameter functions only. These expressions are then fitted 
to all available K — > 2ir and K — > 3ir data to obtain a fit of the parameters that 
occur. We compare with the work of Kambor, Missimer and Wyler. 
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1 Introduction 



Chiral Perturbation Theory (ChPT) is the low-energy effective field theory of the strong 
interactions. It was introduced in its modern form by Weinberg, Gasser and Leutwyler 
|l|, 0, [|. It has had many successes and applications. A pedagogical introduction can 
be found in ||]. The method has been used as well for nonleptonic weak decays. The 
main work of extending ChPT to the nonleptonic weak interaction was done long ago 
by Kambor, Missimer and Wyler who worked out the general formalism || and applied 
it to K — > 27r, 37r decays ||. These results were then used to obtain directly relations 
between physical observables in J7j. Reviews of applications of ChPT to nonleptonic weak 
interactions are ||. 

Earlier work using current algebra methods or tree level Lagrangians relevant for K —>■ 



37r are |], |10[ and references therein. Many new measurements in K — > 3ir have become 
available since the work of |J so an update of the fits of the parameters done there became 
necessary. The analytical expressions of [fj were never published and have since been lost. 
This forced us to reevaluate these decays and we present the analytical results here using 
a simplified form derived using the arguments first used for 7T7r-scattering by Knecht et 
al. ||11|| . We then perform a detailed comparison of the expressions with the data and find 
a reasonable agreement. We point out some directions for future work. 

The next section defines our notation and gives a unique list of possible terms at next- 
to-leading order in ChPT. Section ^ gives the list of decays and the isospin constraints as 
well as the simplified analytical form valid up till order p 6 in ChPT. The ambiguities in this 
simplified form are discussed in detail in App. [A]. This section also discusses the isospin 
constraints on the various amplitudes. Our main result, the recalculation of the K — > 3tt 
amplitudes is discussed in Sect. f| where we present the tree level result and work out the 
independent combinations of the weak p A parameters that occur. The expressions at order 
p 4 are still rather cumbersome and are listed in App. ||. Section |5| contains a short review 
of the available data and presents in detail the various fits we have performed to the data. 
Our conclusions are summarized in the final section. 



2 The ChPT Lagrangian 

The Lagrangian of ChPT in the strong sector was worked out at next-to-leading order 
(NLO) in and is given by 

C s = £ 52 + £s4, (1) 



with 



C S 2 = ^{(D^D^U) + (^U + X ^)}, (2) 

Csa = L^D^D^U) 2 + L 2 (DjfD v U)(D*lf<D v U) 

+Lz(D^DJJD u rfD v U) + L A (D^DJJ)(x ] U + X U ] ) 
+L,{D^D^U( X ] U + rf X )) + L^U + X U ] ) 2 



1 



+L 7 ( X ] U - X U ] ) 2 + L B (xWU + XU ] XU^) 
-iL 9 (F*D»UD u U ] + F^ U D^U ] D U U) 

+L w (tfF*UF L n- 



(3) 



(A) stands for the flavour trace of the matrix A, and Fq is the pion decay constant in the 
chiral limit. The special unitary matrix U contains the Goldstone boson fields 



U = exp {^y-M 



M 
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K~ 



(4) 



7§V 



The formalism we use is the external field method of with s, p, and r M matrix val- 
ued scalar, pseudo-scalar, left-handed and right handed vector external fields respectively. 
These show up in 

X = 2B (s + ip) , 



in the covariant derivative 

and in the field strength tensor 

K» R) = Wu - dj(r), - i P(r)„, Z(r)„] . 
For our purpose it suffices to set 



(5) 
(6) 

(7) 



m d 



p = , 1^ = r M = . 



(8) 



V 



m. 



We also define the matrices u, m m and x±, 

Un = iu'DJJu* =u* u , u 2 = U, x± — u X u =t u\ u . 



(9) 



The ChPT Lagrangian for the weak nonleptonic interactions at lowest order dates back 
to current algebra days M and is 



W2 



CF* 



G 8 (A 32 «X) + G' 8 {A 32X+ ) + G„&> kl (A ijUfl ){A kl u^) 



+ h.c. 



(10) 



The tensor t %h has as nonzero components 



£21,13 _ £.13,21 



£23,33 _ £33,23 



3 
1 



£22,23 _ £23,22 



1 

6 



t 



23,11 _ £H,23 



1 

3 



2 



The matrix Ay is defined as 



A,; 



(A 



V lab 



Oia Ojb 



(12) 



We only quote the |AS| = 1 part here but the others can be obtained by changing the 
indices in A 32 and t^ ,kl appropriately. The parts with G 8 and G 8 transforms as an octet 
under SU(3) l while the G 27 part transforms as a 27 under the same group. The term with 
G 8 is often referred to as the weak mass term. 

The coefficient C is defined such that in the chiral and large N c limits Gg = G27 = 1, 



C 



3 Gp v V* 



The correspondence with the parameters c 2 and c 3 of || § is 

c 2 = CFqG$ c 3 = —-CFqG 2 7 ■ 

o 

The NLO nonleptonic weak Lagrangian was first worked out in 
there was redundant. Subsequent work on a less redundant basis was ]I2|]. A fully nonre- 
dundant basis for the octet part was presented in [13]. We will use the basis of [H| for the 
octet part and use the notation of for the 27 part, but keep only the nonredundant terms 
as derived in |Ll|. The NLO weak ChPT Lagrangian, quoting only the terms relevant for 
K — > 2tt and K — > 3n decays, is 



(13) 
(14) 

but the basis given 



C 



(4) 

A5=l 



C F 2 lG s [NtOl + N 2 Ol + N 3 8 + N 4 0lN 5 8 5 + N 6 8 + N 7 O s 7 

+N 8 O s s + N 9 O s 9 + N 10 O 8 10 + N n 8 n + N 12 8 2 + N 13 8 3 

+G 27 \DiO? + D 2 2 2 7 + D 4 Of + D 5 Of + D 6 Of + D 7 0? + D 2e Ojl 



+D 27 0% + D 28 Of 8 + D29OH + D 30 OH + D 31 OH 
The octet operators are 



+ h.c. 



(15) 



o 8 = 


(A 32 M M u M M i ,M I/ ) , 


o\ = 


(A 3a it M u 1 ,u , V) , 


ol = 


(A 32 u tl u u )(u^u L ') , 


o 8 = 


(A 32 it /1 ) {uyU^u 11 ) , 


o 8 = 




ol = 




o 8 = 


(A32X+>(«/), 


o 8 - 


(A 3a «X><X+) J 


o 8 - 


(A 32 (^ti/-ii/x^)) , 



3 



The 27 operators are 
Of 



Of 

of 
of 

^6 

Of 

o 27 

^26 

Of 

of 

^29 

27 
^30 

^31 



O 8 

O 8 

o 8 

w 12 

o 8 

^13 



t ij,k. 
f ij,k. 
t ij,k. 
t ij,k. 
fj,k. 
t ij,k. 
t ij,k. 
f ij,k. 
t ij,k. 
t ij,k. 
t ij,k. 



= (A32X+X+) , 
= (A 3 2X+)(X+), 
= (A32X-X-) , 

= (A 32 X-)(x->- 

(A ijX -)(AkiX-) , 
(A ij u fM )(A kl {u tl x+ + X+u tl )) , 
(A ijUfl )(A kl (u^ -x-u")), 

(AyX+}(Afci«X) » 

(A ii ^)(A fe ^)( x+ ), 

(AijU^) {A kl u u u u ) , 

(Aij(u^u u + UvUfj)) (A H (uV + uV)) , 

(Aij^Uu - u u u^)){A k i{u^u u - //)) , 

(A ij u fl )(A k iu u u fM u u ) , 

(AijU^) (A k i(u^u v u u + ti/«")> , 

(A ii ^)(A^)(w^). 



(16) 



(17) 



Notice that some have the opposite sign compared to ||, because of the difference in the 
definition of x% an d P of that reference. The basis used in [H| for the octet case is slightly 
different but related to the one here by 



N 5 ~- 
N 36 



E10 — E n 

En ; 
~-E 2 - 
'-E h . 



N 9 -- 
N 12 



En + 2E12 
E15 ; 

= -E 3 + E 5 : 



N 7 = 
N 13 



I -En 
-E 1 - 
- Ei ; 



'r E 13 ; 
E5 ; 



Notice that there is a small misprint in the relation between the Ei and iVj in (2.19) of |14[ 
in the relation between Ni 3 and E4. 

The infinities appearing in the loop diagrams are canceled by replacing the coefficients 
in ([H) and (|T5D by the renormalized coefficients and a subtraction part. The infinities 
needed in the strong sector were calculated first in |J and those for the weak sector in || 
and confirmed in [0. For the terms in Eqs. (§) and fll~5l ) they are all of the type 

1 



X: 



(e» 



-2c 



16n 2 e 



with the dimension of space-time d = 4 — 2e and 



c=--(ln(47r) + r'(l) 



(19) 



(20) 



for X = L,N, D. The coefficients are listed in Table [I] || |5|, |T2" ]. 
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N- 




Tli 




1 
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3 /fi4 
— o/ U4 


i 
i 




1 

— 1 
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1 /1 9 
1/ 1Z 


o 
L 


— o/ OZ 


o 
L 




1 /zL 
1/4 


o 
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n 
U 


Q 
O 


n 

U 


Q 
O 




n 

U 


A 
4 


— o/z 


A 
4 


1 /I P. 
— 1/ 10 


A 
4 




1 /9 
-1/2 





1 /9 
-1/Z 


K 
o 


— ^ /1 (S 
O / 1U 


c; 


— ^ /4 - 

9/4 


_ Q/Q f/^/ /(7 "1 


u 


^/4 
9/4 


6 


—1 1 /?88 


6 




1/8 


7 


-1/2 


7 





7 


9/16 - 


- 1/4 (G' 8 /Gs) 


26 


1/2 


8 


-5/96 


8 




1/4 


27 


1/4 


9 


-1/8 


9 


-3/8- 


h 3/8 (G' 8 /G 8 ) 


28 


5/6 


10 


1/8 


10 


-1/3- 


5/24 (G' 8 /G 8 ) 


29 


-19/ 






11 


13/36 - 


11/36 (G 8 /G 8 ) 


30 


-5/3 






12 


5/24- 


5/24 (G 8 /G 8 ) 


31 
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Table 1: The coefficients of the subtraction of the infinite parts defined in Eq. fll 

3 Kinematics and Isospin 

There are five CP-conserving decays of the type K 



3tt: 



K L {k) - 




[^000 ] J 


K L (k) - 






K s (k) - 


-> 7r + (pi)7r"(p 2 )vr (p 3 ) , 




K + (k) - 


- ^(Pi)n°(p 2 )^(p 3 ), 


[Aoo+] 


K + (k) - 


-> 7r + (pi)7r + (p 2 )vr _ (p 3 ) 


[A ++ - 



(21) 

where we have indicated the four momentum defined for each particle and the symbol we 
will use for the amplitude. The K + decays have an obvious counterpart in K~ decays. 
The kinematics is normally treated using 

s 1 = (k-p 1 ) 2 , s 2 = (k-p 2 ) 2 , s 3 = (k-p 3 ) 2 . (22) 

The amplitudes are often expanded in terms of the Dalitz plot variables x and y defined 

S 2 — Sl S3 — Sq 1 / 2 2 2 2 \ , . 

S = o \ m K + "V + "V + "Vj » ( 23 ) 



via 



.r 



2/ 



2 

m + 



where the kaon mass and the pion masses in sq are those from the particles appearing in 
the decay under consideration. 

The amplitudes for Ki — > 7r + 7r~7r°, K + — > 7r + 7T + 7r~ and — > 7r°7r°7r + are symmetric 
under the interchange of the first two pions because of CP or Bose-symmetry. The ampli- 
tude for Kl — > 7r°7T 7r is obviously symmetric under the interchange of all three final state 



5 



particles and the one for Kg — > 7r + 7r 7r° is antisymmetric under the interchange of 7r + and 
tt~ because of CP. 



Isospin invariance does give some constraints on the amplitudes [|15], [16], [D| . These can 
be written as 



^ooo ( s i > s 2, S3) 

A+_ (si,S 2 ,S 3 ) 
A^_ (si,S 2 ,S 3 ) 
Aoo+(S!, S 2 , S3) 

A ++ -(si, s 2 , s 3 ) 



3A [Si, S2, S3) , 

A n (s h s 2 , s 3 ) - B n (s 1: s 2 , s 3 ) , 
2 

Co(si, s 2 , s 3 ) + - [-B t (s 3 , s 2 , si) - B t (s 3 , s h s 2 )} 

A c (si, s 2 , s 3 ) - B c (s!, s 2 , s 3 ) + B t (st, s 2 , s 3 ) , 
2A c (s!, s 2 , s 3 ) + B c (si, s 2 , s 3 ) + B t (st, s 2 , s 3 ) . 



(24) 



The functions A c ^ n are fully symmetric in si, s 2 , S3. Co is fully antisymmetric and B c>n j(si, s 2 , S3) 
is symmetric under the interchange si <-» s 2 and satisfies the relation 



5j(si, s 2 , s 3 ) + 5i(s 2 , s 3 , si) + Bi(s 3 , si, s 2 ) = . 



(25) 



74 cn and B cn belong to the 1 = 1 final state, B t to the 1 = 2 final state and Co to the 
1 = final state. 

A further simplification can be obtained by observing that the imaginary parts of loop 
diagrams in ChPT do not contain I > 2 until at least 0(p 8 ). This allows to rewrite the 
amplitudes in a way that only contains functions of single variables and is fully correct up to 
0(p 6 ). The underlying arguments are the same as the ones used for a similar decomposition 
for 7T7T scattering by ||11|| . 

A£ 00 (si,s 2 ,s 3 ) 
A L 



ySl, S 2 , S 3 J 



Aj_ (si,s 2 ,s 3 ) 

A 00+ (si, s 2 , s 3 ) 
A ++ „(si,s 2 ,s 3 ) 



M (si) + M (s 2 ) + M (s 3 ), 

Mx(s 3 ) + M 2 ( Sl ) + M 2 (s 2 ) + M 3 (s!)(s 2 - s 3 ) + M 3 (s 2 )(s! - s 3 ) , 
M 4 (si) - M 4 (s 2 ) + M 5 (si)(s 2 - s 3 ) - M 5 (s 2 )(si - s 3 ) 
+M 6 {s 3 )(s 1 - s 2 ) , 

M 7 (s 3 ) + M 8 (si) + M 8 (s 2 ) + M 9 (s!)(s 2 - s 3 ) + M 9 (s 2 )(si - s 3 ) , 
M 10 (s 3 ) + M u ( ai ) + M n (s 2 ) + M 12 (sx)(s 2 - s 3 ) + M 12 (s 2 )( Sl - s 3 ) . 

(26) 



The expressions for the various amplitudes become significantly shorter when written ex- 
plicitly in this form. Comparing (p4[) and ([26]) gives 

A>(si,s 2 ,s 3 ) = i M o(si), 



i=l,3 



^4 c (si,s 2 ,s 3 ) 

B n (si, s 2 , S3) 
5 n i(s) 



£ - (M 7 (s 4 ) + 2M 8 (s,)) , 



i=l,3 



S„l(si) + ^nl(s 2 ) - 2 J B„i(s 3 ) + S n2 (s 1 )(s 2 - S3) + £ n2 (s 2 )(Si - S3) 

±( Ml (s)-M 2 (s)) , 



B n2 (s) = -M 3 (s), 

B c {sx, s 2 , s 3 ) = S cl (s x ) + S cl (s 2 ) - 2-B cl (s 3 ) + S c2 (si)(s 2 - s 3 ) + 5 c2 (s 2 )(si - s 3 ) , 

fi c i(s) = ^(M 7 ( S )-M 8 ( S )-M 10 ( S ) + M n ( S )) 
o 

5 c2 (s) = i(-Af 9 (a) + Mi 2 (s)) , 

5 t (si, s 2 , 8 3 ) = B tl ( Sl ) + B tl (s 2 ) - 2B tl (s 3 ) + B t2 ( Sl )(s 2 - s 3 ) + 5t 2 (s 2 )(si - s 3 ) , 

B tl (s) = ^(-M 7 ( S ) + M 8 ( S )-M 10 (s) + M n ( S )) 


S t2 (s) = ^(M 9 (s) + M 12 (s)) , 
C (si,s 2 ,s 3 ) = Coi(s 3 )(si - s 2 ) + C 01 (si)(s 2 - s 3 ) + C i(s 2 )(s 3 - s x ) , 

Coi(s) = i(2M 5 ( S ) + M 6 ( S )) (27) 

and the relations 

M (s) = M 1 {s) + 2M 2 {s), 
2M 7 (s) +4M 8 (s) = M 10 (s) + 2M n (s), 

M 4 (s) = ^(M 7 (s)-M 8 (s) + M 10 (s)-M n (s)) 
M 5 (s)-M 6 (s) = M 9 (s) + M 12 (s). (28) 

The split-up between the polynomial parts in (ffiS|) has some ambiguity as discussed in 
App. The same ambiguity appears in the rewriting of the B cnt functions in single 
variable parts and can affect the last two relations in fl2"E|). The results given explicitly in 
App. [B| are brought in a form satisfying the relations ( p8[) using this freedom. 



4 Analytical Results 

The lowest order result is well known and follows from the diagrams in Fig.|l|. The Feynman 
amplitude at tree level for Aq 00 is 

C 1 

M o(^ = ^p^^(G 8 -G 27 )ml. (29) 

A choice for A+_ compatible with the relations (pS|) and (|29| ) is 

CP 4 1 

M 1 (s)\ p2 = ijL±-(G B -G„)m 3 K , 
M 2 (s)\ p2 = 0,* 

Ms(s)\, = .•0-{4 + ^=5^4) (30) 
p F*F K I 3 6 - mi 
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Figure 1: The tree level diagrams for K —>■ 3tt. A filled square is a vertex from Cw2 and a 
filled circle from C$2- 

The tree level result for A+_ allows very much of reshuffling between M 4 , M 5 and M 6 . A 
choice is 



s 



<p F£F K mr K - mi 6 v 

Me(s)\ p2 = i^^^l^ml-Sml). (31) 



A 00+ leads to 



A' 



G*27 5_ 2 _2 , 13_ 4 1 _4 , 7_ 2 17 



mi — ml V 3 o 2 o o 



M 8 (,)| p2 = -M 7 ( a )|^+»-^-f- 2 ^ 27 - £ (2m^ - 3m^) s , 
M 9 (s)\ p2 = 0. (32) 
The amplitude for can be written as 

M w (s)\ p2 = i^j- {G 8 (s-m2-m^)+Gw(-yS + m|- + ym£)} , 
M n (s)| p2 = 0," 

M 12 (s)L = 0. (33) 



The result at order p 4 is significantly longer. The diagrams that contribute are shown in 
Fig. |2|. Using the decomposition of all the amplitudes in the Mj(s), and using the freedom 
in choosing the functions Mj(s), relatively compact expressions can be obtained and we 
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D— O 



■ 6 




Figure 2: The diagrams of order p 4 . An open square is a vertex from Cwa, an open circle 
a vertex from £54, a filled square a vertex from Cw2 and a filled circle a vertex from £$2- 
Diagrams contributing to strong interaction wave function renormalization are not shown. 
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T." 


S~1 f ATT O ATT 1 O ATT 1 ATT\ 1 /T f 1 7"Yr \ 

G 8 (A 5 - 2A 7 + 2A 8 + N 9 ) + G 27 {~ 2 D 6 ) 


T ' 

K 2 


G 8 (N{ + N r 2 ) + G 27 (±D r 26 - J 


K 3 


G 8 {ND + G 27 [\D r 27 + lD r 28 ) 


K± 


G 27 {pi - Dl + 


K 5 


G 27 (D^o + -D31 + 2D$g) 


K G 


G 27 (81% - £% + ^0) 


K 7 


G 27 (-41% + -D29) 


Ks 


G 8 (2A£ + 4iVf + jV£ - 2iV[ - 4iVT a - 2iV{ 2 ) + G 27 (-§£>[ + fDg) 


K9 


G 8 (iV 5 r + ^ + ^) + G 27 (-| J Dg) 


Kio 


G 27 - 2L>£ - D?) 


Ku 


G27 



Table 2: The independent linear combinations of the N[ and Z?J that appear in the am- 
plitudes. 

have given them explicitly in App. |B[ Our results for K —>■ tttt are in full agreement with 

rj . 

In order to do the data fits it is important to see how many combinations of the 25 
unknown parameters appearing in fllBT) actually matter. Taking the various independent 
coefficients from the amplitudes Mj(s), it turns out that only 11 linear combinations of the 
25 parameters show up. Of these seven appear in the amplitudes multiplied by coefficients 
of order m A K while the remaining four have coefficients of order m 2 m 2 K . 

The measurable combinations are displayed in table |2|. Any linear combination of 
them is of course also allowed. The choice displayed in the table was made to have all 
combinations to be orthogonal to the non-measurable ones when possible. The remainder 
was then chosen simply to have somewhat compact expressions for their contributions to 
the Mj(s) and the A Q and A 2 amplitude for K — > txtx. K4 is the combination multiplying 
m A K in the A 2 amplitude. K\ is the octet part of the A amplitude proportional to m A K plus 
the part of the 27 multiplying m A K that cannot be written in terms of K4. K\q and K§ are 
the equivalent parts of the coefficients of m\ in A 2 and Aq. The numbering is chosen such 
that K\ to K 7 contribute leading in m 2 K and inside those categories the ones with octet 
contributions are listed first. There are 5 combinations with octet contributions of which 
three leading in m 2 K . 

The amplitudes we calculated are at one loop in ChPT. They thus include final state 
rescattering effects. A partial analysis of this type of effects was done in ]I7J. They derived 
the relations 

/ A c (s 1 ,s 2 ,s 3 ) \ _ , . ( A c (s u s 2 ,s 3 ) \ 

{B c ( Sl ,s 2 ,s 3 ) ) R ~ R{SuS2 > S3) {b c ( Si ,s 2 ,s 3 ) ) NR 

\ B n ( Sl ,s 2 ,s 3 ) J R K x ' 2 ' 3J \ B n { Sl ,s 2 ,s 3 ) ) NR 
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Figure 3: Im(v4 ++ _) as a function of x and y normalized to — |A ++ _(s , s , sq)\. The 
curve with the plateau at zero indicates the phasespace. The thresholds in the various 7T7T 
channels are clearly visible. 



B 2 (s 1 ,s 2 ,s 3 ) R = 5(s 1 ,s 2 ,s 3 )B 2 (s 1 ,s 2 ,s 3 ) NR 



(34) 



for the amplitudes defined in (0). The subscript R on the l.h.s. means that rescattering 
effects have been included and NR on the r.h.s. that they are not included. The expressions 
for the matrix function R(s\, s 2 , S3) and the function S(si, s 2 , S3) to lowest order in ChPT 
can be found in |I7| |. We have checked numerically that our amplitudes satisfy those 
relations to the order required. The approach to the various thresholds (sj — > 4m^) is as 
expected to this order in ChPT. The effects of the thresholds can be seen in Fig. [3[ The 
plot shows Im(A ++ _) as a function of x and y normalized to — \A ++ _(s , s , s )\ 



5 Experimental Data and Fits 



At the end of the seventies a review was written incorporating the then finished precision 
experiments on the charged kaon decays ||18|| . Since then there have been many new results 
on K — > 37r decays. A recent attempt to refit the data using some of the partial numerics 
of is fl"9|l . We have in this paper recalculated the full expressions and will do a full fit 

|T5| . We analyze the data in terms of the widths given in Table [5]. 

We do not directly use 



to the data contrary to 



20 



The numbers are taken from the review of particle properties 
the data for r^_ but instead fit directly to the measured quantity A as discussed below. 
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Decay 


Width [GeV] 


CHPT fit [GeV] 


K S -> 7T 7T 
— > 7T + 7T~ 


(1 1 245 ± f) 0078") • 1 -17 
(2.3124 ±0.0207) • 10~ 15 
(5.0543 ±0.0211) • 10~ 15 


1.124 • 10~ 17 
2.312 • KT 15 
5.054- KT 15 


K L -> 7r 7T 7T U 

if/, — > 7r + 7r~7r° 
i£s — > 7r + 7r~7r° 

K + -> 7T 7r 7T + 
— > 7T + 7r + 7r~ 


(2.6901 ±0.0402) • 10~ 18 
(1.5978 ±0.0283) • 10~ 18 

(2.36 ± 0.81) • 10- 21 
(0.9194 ±0.0426) • 10~ 18 
(2.9706 ±0.0272) ■ 10" 18 


2.655 • KT 18 
1.626 • 10~ 18 
3.1 • 10- 21 
0.911 • 10~ 18 
2.973 • 10~ 18 



Table 3: The various decay widths from the PDG tables|20| . 



In addition several slope parameters are measured in the various decays. The distri- 
butions in the Dalitz plot are conventionally described in terms of x and y defined in 



The decay amplitudes squared are now expanded as 



A(si,s 2 ,s 3 ) 



i + gy + hy 2 + kx* 



(35) 



where we have used CP invariance and the symmetries in the decays. For Kl — > 7r°7r 7r° 
g = and k = h/3. For Ks — > 7r + 7r~7r° one defines the ratio A or measures directly the 
expansion of the amplitude ET| 



Umax f-^lim 

(f) 

dy / dxA^_ A s + _ 

Vmin ^0 

Umax fXlim(y) . T 

dy / dxA^*_ Al_ 

Vm in ** 



75a; - isxy . 



(36) 



The various measurements are given in Table [|. Experiments after 1980 have been cited 
explicitly. Earlier references can be found in the PDG[^(| or in the comprehensive re- 
view ||18|| . Notice that there are significant discrepancies between the various experiments 
indicated by the fairly large scale factors. 

The data have been analyzed previously in terms of expansions in x and y of the 
amplitudes. 



4oo = 3(a 1 ±a 3 )±3(Ci-2C 3 )(l/ 2 ±^ 2 ) 



^+-0 = («i + «s) - (ft + fo)y + (& - 2( 3 ) [y 2 + ~x 2 j + (£1 - 2£ 3 ) [y 2 - -r 

2 4 
^+-0 = gV373»- g^ 2 ^) 
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Decay 


Quantity 




s 


Ref. 


CHPT fit 


Kl ^ ^u^u 


ft, 
ft 
h 


-0.0033 ±0.0013 
-0.0061 ±0.0010 
—0 00506 ± 00135 


1 7 


I 

1 

av< 


22 


ige 


—0 0072 


Kl — > 7r + 7r _ 7r° 


9 
h 

k 


0.678 ± 0.008 
0.076 ±0.006 
0099 ± 001 5 


1.5 




21 
21 
21 


• - 


m 


0.677 
0.085 
0055 


Ks — > 7r + 7r _ 7r° 


Re(A) 
Im(A) 


0.0316 ±0.0062 
-0.0088 ±0.0068 
f3 3 + Ft) ■ 1 0~ 8 
(0.4 ±0.8) ■ 10" 8 






|2l 

21 

f 


• - 

• • 

21 
21 


24} 


0.0359 
-0.003 
S 4 • 1 0~ 8 
-0.2 ■ 10~ 8 


K ± _^ 7r u 7r u 7r ± 


9 
h 
k 


0.652 ±0.031 
057 ± 018 
0.0197 ±0.0054 


2.7 
1.4 


w 

2, 


% 

). 

21 


26 
(' 


, 2C] 
2J 


0.638 
0.074 
0.0045 


K + — ► 7T + 7r + 7r _ 


9 
h 
k 


-0.2154 ±0.0035 

0.012 ±0.008 
-0.0101 ±0.0034 


1.4 
1.4 
2.1 


1 


F° 

20 
20 


I 


-0.216 
0.012 
-0.0052 




9 
h 

k 


-0.217 ±0.007 
0.010 ±0.006 
-0.0084 ±0.0019 


2.5 
2.1 


I 
1 
I 


20 
20 
20 


1 





Table 4: The measurements of the Dalitz plot distributions. Notice that in many cases 
the scale factors calculated with the PDGj20| method are sizable. The references refer 
to the most recent measurements and the PDG. S is the scale factor calculated using the 
prescription of PU . 



A)o+ = (-«i + -as) + (pi - ^ 3 - v^ts) y + (-G - Cs) (y 2 + 

A++- = {-2a, + a 3 ) + (-fa + ^ 3 - v 7 ^) V + (~2Ci - 2(a) (y 2 + \x 2 ^j 

+{Hi + ^-Q(y 2 -\x 2 ) • (37) 

A main problem in dealing with the data fitting is the question of how to deal with the 
isospin breaking effects. The phase-space itself is very sensitive to the precise masses of the 
pions and kaons which are used. There are several ways to deal with this problem. One is 
to first fit the expressions ( |37"D to the data where x and y are defined with so calculated 
with the physical masses but otherwise the charged pion mass as in (|23|) . The result of 
this fit is shown in Table [5[ Here we have assumed that the only complex phase appearing 
is the relative phase between A and A 2 . The errors quoted are the minuit errors. For the 
i^-measurements only Re(A) is included. Leaving it out changes the y^/DOF to 3.9/4 
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quantity 


Ref. [ 


18 1 


Ref. | 


1 


^jur nx 


p2 


Aj — U 


1 A 1 

\ A o\ 


0.4687 ± 0.0006 


0.4699 ± 0.0012 


U.40ZZ ± U.UU14 


input 


input 


1 A 1 
\ A 2 1 


0.0210 ±0.0001 


r\ nm 1 i r\ r\r\r\-\ 

0.0211 ± 0.0001 


U.Uzlz ± U.UUUl 


input 


input 


X X 
2 — 


(-45.6 ±5)° 


(—61.5 ± 4)° 


( CO O _J_ /I \0 

— Do.Z ± 4 J 








91.4 ±0.24 


91.71 ± 0.32 


yo. ID ± 0.60 


1A C\('7 r l K\ 

( 4.U( ( 6.0 ) 


KO A 

oy.4 


a 3 


-7.14 ±0.36 


— ( .60 ± U.4 < 


79 I n A fi 
— 0. ( Z ± U.4D 


A Q( A Q\ 
— 4.0^4.0 ) 


— 0.0 


Pi 


-25.83 ±0.41 


OK «Q _(_ n 97 


_97 rifi 4- n 


— 1 7 7(1 (\ 9\ 


—91 Q 


H3 


-2.48 ± 0.48 


-2.43 ±0.41 


— 2 22 ± 47 


— 1 2(1 1) 


— 1 


93 


2.51 ±0.36 


2.26 ±0.23 


2.95 ±0.32 


2.3(2.1) 


2.5 


Ci 


-0.37±0.11 


-0.47 ±0.15 


-0.40 ±0.19 




0.26 


Cs 






-0.21 ±0.08 


-0.09 ±0.10 




-0.01 




-1.25 ±0.12 


-1.51 ±0.30 


-1.83 ±0.30 




-0.46 








-0.12 ±0.17 


-0.17±0.16 




-0.01 








-0.21 ±0.51 


-0.56 ±0.42 




-0.06 




12.8/3 


10.3/2 


5.4/5 







Table 5: The values of the parameters in the expansion in x and y of the K — > 6n 
amplitudes and the K — > Tin amplitudes. The \ 2 quoted are for the fits with the then used 
data. \Aq\ and \A%\ are in units of 10~ 6 GeV . ai,. . . ,£3 are in units of 10~ 8 . 

but all parameters stay within the errors given. 

In the column labeled p 2 in Table |5] we fixed Gg and G27 from the result for A Q and A 2 
at tree level and show the predictions for the linear slopes in the amplitude at this order. 
The source of the differences with [Q is not obvious but could be due to a different pion 
and kaon mass. The numbers in the table were calculated using the charged pion and kaon 
masses, using the neutral ones instead leads to the numbers in brackets. As can be seen 
a general reasonable agreement within about 40% is obtained. The input here determines 
G$ = 10.4 and G27 = 0.61 when we used for F the pion decay constant F n and used F v 
and Fk in the tree level amplitude with the physical values. 

Going to the next order, we first show also in Table [| in the column labeled Ki = the 
p 4 results with G$ and G27 fixed to fit \A \ and \A 2 \. These were calculated with charged 
pion and kaon mass. The input here determines G§ — 5.47 and G27 = 0.39. The change 
compared to the previous is an indication of the size of the p 4 effect. It should be kept 
in mind that we have chosen to normalize the lowest order including one factor of 1/F^, 
changing that to 1/F n changes the relative size of the p 2 and p 4 contributions significantly. 
There is as well a mild dependence on the value of the L[ used. We use in general a scale 
of fi = 0.77 GeV and the L\ values of the one-loop fit of f27], the p 4 fit including the latest 
Kg4 data. These values are (all at /i = 0.77 GeV.) 

L\ = 0.38 ■ 10~ 3 , L r 2 = 1.59 ■ 10~ 3 , L r 3 = -2.91 • 10~ 3 , L\ = 0, 

L\ = 1.46 ■ 10~ 3 , Lg = 0, Lj = -0.49 ■ 10" 3 , L r 8 = 1.0 ■ 10~ 3 . (38) 

It is not possible to determine all 13 quantities G$, G27 and K lr . . ,K n . In principle it 
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cannot be done if we try to fit to the 12 quantities listed in Table |5|,[] but in practice the 
situation is even worse, simply fitting the possible parameters without constraints leads to 
fits with negligible tree level contributions. We have therefore chosen a two-step process, 
first we study the variation of fitting subsets of G$, G27 and K\ r . . ,Kn to \A \,. . . ,£ 3 . 
This allows for a more direct study of the various dependencies on parameters. Finally we 
present the direct fit where we compare directly to the data, but allowing in addition for 
a variable phase 62 — Sq between A and A 2 

The various relations advocated in [[F] follow if one sets 

K x = K A = and K 8 = ■■■ = K u =0. (39) 

The resulting fit values are in Table || column 2. Several of the Ki are uncomfortably large 
but inspection of the size of the p 2 versus p 4 contributions shows nothing conspicuous. The 
source of the problem turned out to be 73 which with the constraints of Eq. (|39"D can only 
be fitted by varying K e but it comes multiplied with a small factor. It then produces fairly 
large other Ki in order to minimize the effect of K 7 on the other quantities. We have 
therefore also fitted using the constraints 

K x = K 4 = K 6 = and K 8 = ■ ■ ■ = K u = . (40) 

We can now attempt to also determine somewhat the combinations K\ and K4. Leaving 
Ki free hardly changes the quality of the fit but shows that the actual value of G 8 can 
change over a fairly wide range. The results are shown in the fourth column of Table |6| 
with constraints 

K 4 = and K 8 = ■ ■ ■ = K n = . (41) 

The same type of correlation occurs for G27 and K4 as shown with the results in column 
5 of the Table where we fitted with constraints 

K 1 = 0, K 4 = 0.1 and K 8 = ■ ■ ■ = K u = . (42) 

The reason we have picked a fixed value for K 4 is because there is very long, narrow and 
extremely shallow fitting region that moves in the end to a very small G27 and values for 
K4, K5, Kq and Kj that are simply enormous but a total y 2 which is only marginally 
smaller than the one shown, 11.7 rather than 11.9. 

How well do the relations proposed in |7| stand up to the new data. In practice, for the 
octet which dominates the AI = 1/2 part, it means that oti determines (1, (3\ determines 
£1 . The first of these relations is extremely well satisfied, but it is hard to get £1 to better 
than about 2.5 standard deviations, the value obtained from the various above fits is 

6 w (-0.98 -i0.13) • KT 8 . (43) 

Inspection of the numerical coefficients for the octet parameters that only contribute sup- 
pressed by powers of m 2 reveals that some of them show up with large numerical coeffi- 
cients. One can then instead choose to leave K 8 or Kg free rather than K\. The resulting 

1 &2 — <$o cannot be fitted at this order. 
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constraints 


Eq. (£9j) 


Eq. ([40D 


Eq. (HU) 


Eq. ([42D 


Eq. (H) 


Fitted 


\Ao\,---,& 


experiment 


Gg 


5.47(2) 


5.47(2) 


7.24(2.04) 


5.49 


5.45(2) 


G 27 


0.392(2) 


0.392(2) 


0.392(2) 


0.139 


0.392(2) 


lO^x/Gs 


= 


= 


-8.5(7.5) 


= 


= 


10 3 K 2 /G 8 


54.7(2.8) 


53.6(2.7) 


41.3(11.4) 


53.7 


51.9(3.2) 


10 3 K 3 /G 8 


3.0(1.4) 


3.5(1.3) 


10.0(6.3) 


3.2 


3.8(1.5) 


10 3 K 4 /G 27 


= 


= 


= 


= 100 


= 


10 3 K 5 /G 27 


-54.5(23.4) 


-19.8(9.6) 


-54.5(23.1) 


-49.4 


-42.5(16.6) 


W 3 K 6 /G 27 


-185(114) 


= 


-184(113) 


-366 


-166(113) 


10 3 K 7 /G 27 


114(46) 


45.1(18.2) 


114(46) 


199 


120(32) 


X 2 /DOF 


12.3/5 


14.9/6 


11.8/4 


11.9/4 


26.8/10 



Table 6: The results for Gg, G 27 and the K% for the various constraints discussed in the 
text. In brackets are the MINUIT errors. 

fits are of similar quality to the one with K\ free and do not significantly improve the 
prediction for £1 given in ( fj3|) . The possible variation is less than 10%. 

There is no similar single value with significant discrepancies for the AI = 3/2 parame- 
ters. The above named problem with 73 does not lead to a major discrepancy but typically 
the quadratic quantities are the ones with a deviation of about 0.5 to 1.4 a. 

Varying the masses of the pion and kaon or the eta mass used in the loops does not 
change the results of the fits significantly. 

The fit via the intermediate step of \A$\, . . . , £ 3 is fairly fast and allows an easy study 
of the variation with the inputs. We have also done the fit directly to all the experimental 
data listed in Tables and |] with only the use of Re(A) for the K$ — > /T + 7r~7r° data. The 
masses used in the phase space were the physical masses occurring in the decays. The 
masses used in the amplitudes are the physical kaon mass and the pion mass such that 
3m% = J2i=i,3 m 1i is satisfied. The eta mass in the loops we then calculated using the 
GMO relation. The resulting values for the are given in the last column of Table [3] and 
the values for the widths and Dalitz plot distribution variables are in the column labeled 
CHPT fit in Tables |3] and |j. We have allowed an extra phase between Aq and A 2 in this fit. 
Otherwise this fit corresponds to the constraints ( j39|) . The higher y 2 is mainly due to the 
fact that £x has been fitted to more than one quantity, and here we have the discrepancy 
with the slope k in — > tt + 7t~7t° and K + — ► 7r°7r°7r + . These together account for 16.4 of 

X 2 - 

The fact that the last fit and the fit via the intermediate step of \Aq\, . . . , £ 3 agree very 
well is an indication that the isospin breaking effects at the amplitude level are small. 
An estimate can be done by calculating the amplitudes and the slopes with the masses of 
K°, 7T° and comparing them with the amplitudes calculated with mjf+, m 7r +. Since often 
factors of m 2 K — appear this is the worst case. The changes are typically 2% for the 
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amplitudes, 5-8% for the linear slopes and up to 20% for the quadratic slopes. 



6 Conclusions 

In this paper we have recalculated the K — > tttt and K — > tttttt amplitudes to next-to- 
leading order in CHPT in the isospin limit and presented the amplitudes in a simple form 
suitable for dispersive estimates of higher orders. 

We have performed a new global fit to all available kaon data at this level taking into 
account the isospin breaking in the phase space. Work is under progress to include isospin 
breaking also in the amplitudes. 

At present the situation is that a satisfactory agreement is obtained with the predictions 
keeping in mind that for the quadratic slopes higher order corrections could be sizable as 
happened with the analogous results in r\ — > 3n results. A discussion of the latter together 



with references can be found in ||28|| . Some work exists on dispersive corrections in K — > 37r 
17| , but more study is needed to fit their results in this framework. 

Now the isospin breaking corrections, including possible electromagnetic radiative cor- 
rections, need to be studied to see if they are the cause of the differences in the quadratic 
Dalitz plot parameters. Work exists to a large extent for the K — > 2tt decays p9|| , but 
little has been done using modern methods for K — > 3n decays. 
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A Ambiguity in the definition of Mi 

The functions Mj(s) defined in Eq. (fffi) are not unique for two reasons. The variables 
si,S2 and S3 are not independent, they satisfy the relation (f23|), and low powers of si,s% 
and S3 can be fitted in more than one of the Mj functions. 
As a simple example, we could have added 

SM (s) = a(s-s ) , (A.l) 

with a any constant, without changing the result, for Aq 00 . The full set of these ambigu- 
ities can be derived by checking the total number of independent terms that exist in the 
polynomial expansion of the amplitudes and comparing it with the the same expansion of 
the Mi(s). 

The changes of the functions that leave all amplitudes unchanged are 
M (s) — > M (s)+6 1 (s-s ) (A.2) 
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Mil 

M x i 
M x i 
Mil 
M 3 , 

Mil 

M 3 , 
M 4 i 
M 4 i 
M 5 , 
M 4 i 
M 5 , 
M 4 , 
M 6 , 
M 4 i 
M 6 , 



The triplets 
leading to a 



Mi 

Mi 
Mi 
Mi 
M 3 

Mi 

M 3 
M 4 
M 4 
M 5 
M 4 
M 5 
M 4 
M 6 
M 4 

M 6 



+ <y 2 , m 2 (s) 

+ 5 3 (s - s )) 
+ 25 4 s, M 2 (s) 

- 2<5 5 s (s - 3s ) . 
+ 5 5 s; 

+ {s - 3s ) 3 

+ S 6 s 2 ; 
+ 5 7 ; 

+ 5 8 s, M 5 (s)- 

+ Sg, M 6 (S) — 

- 5 10 s (s - 3s ) . 
+ 5 n s, M 6 (s) - 



->M 2 (s)- 
M 2 (s) — 
— > M 2 (s) 
, M 2 (s)- 

M 2 (s) 



M 2 (s) - ^5 6 s 2 (s - 9s ) 



-> M 5 (s) + 5 8 ; 
M 5 (s) + 5 9 ; 
M 5 (s)^M 5 (s 
-> M 5 (s) + (Jus; 



+ 5 12 (s 3 - 9s s 2 + 18s 2 s) , M 5 (s 
- 25i 2 s 2 ; 
+ S 13 3s (i 



9s s 2 + 18sqS 



25 13 (s 2 - 9s s 2 



\h ; (A.3) 

>M 2 (s) + 5 3 ( S - So )) ; (A.4) 
- 5 4 s, M 3 (s) — > M 3 (s) + 5 4 ; (A.5) 
^M 2 (s) + 5 5 s( s -3s ), 

(A.6) 

(A.7) 
(A.8) 
(A.9) 
(A.10) 
(All) 
(A.12) 

(A.13) 
(A.14) 



♦ M 5 (s) + 5 12 s 2 ; 
M 5 (s) — > M 5 (s) + <W 



(M 7 ,M 8 ,M 9 ) and (M 10 , M n , M 12 ) have an ambiguity similar to (|A~^)-(|A~^), 
total ambiguity of 23 free parameters. 



B Explicit expressions for the M$(s) in K — ► 37T. 

We write the Mj-(s) functions defined in Eq. (EBf) as 



M(s) = M l (s)\ p2+ ^^f \^ 2 Mf(s) + ||Mf (s) + ^M 2 \s)\ . (B.l) 

The effect of G' s cannot be distinguished from higher order coefficients in decays not 
involving external fields. This was known at tree level earlier and has been proven to 
one- loop in ||. The effect of Mf'(s) can be reconstructed from Mf (s) by the changes 

A[ + 2N r u + A[ 2 — A[ + 2JVTi + ^2 + (2^5 " ^ - 8L r 8 ) , 

(jrg 

N r 7 — > JVf + S(_4LJ). (B.2) 

The results for the X — > 7r7r amplitudes are in full agreement with [TJ]] and can be 
found there. 
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We have brought the M/(s) in a form that satisfies the isospin relations (|28|). These 
can be used in the form 

M{{s) = Ml{s)-2Mi{s), 

M>(s) = Mi(s) - M£(s) - M( 2 (s) , 

Mi{s) = 2Mi(s)-±M{ (s) + Mi 1 (s), 

Ml{s) = -M{{s) + ~M{ (s) (B.3) 

to reconstruct these. We have extensively used the GMO relation in writing them in this 
form. 

The explicit expressions for A, B and Bi, the finite part of the loop functions, can be 



found in many places, e.g. |30 
The octet ones are: 



+ (iV[ + iV 2 r ) (-2mlm 2 K + 3ml s - 2ml + m 2 K s - s 2 ) 

+^"3 (2 ml m\ - 3m^s + 5m* - m 2 K s + m A K — 2 s 2 ) 

+N; (2 ml m\ - 2 ml s + 2 m\ + 2 m 2 K s) + N r 7 (8 ml m 2 K - 4 m\ s) 

+N; (-2mls + 2ml + Am 2 K s)+ iV 9 r (-2 ml m\ - 2 ml s + 2 m\ + 2 m\ s) 

+ {N[ Q + 2N[ l + N[ 2 ){-Amlm 2 K ) 

+L\ (32 ml m 2 K -48mls + 32 m% - 16 m\ s + 16 s 2 ) 

+L r 2 (-16 mlm 2 K + 24m 2 n s - 40ml + 8m 2 K s-8m 4 K + 16 s 2 ) 

+L T 3 (16 mlm 2 K - 24m 2 n s + 16ml - 8m 2 K s + 8 « 2 ) 

+L4 (—16 ml m 2 K + 16 ml s — 16 ml) + L\ (—8 ml m 2 K + 16 ml s — 16 m n ) 
+ (l/(16vr 2 )) (4/3 ml m 2 K + 1/2 m 2 n s - 1/2 ml - 3/2 m 2 K s + 1/2 m%) 
+A(ml) (-19/8 ml - 2/3 m 2 K + 21/8 s) + A(m 2 K ) (-13/4 ml - 1/3 m\ + 3 s) 
+A(m 2 v ) (-3/8 ml + 3/8 s) + 5(m 2 , m 2 , s) (1/2 m>^-2m^ + mt + s 2 ) 
+5(m 2 , m 2 K , s) (13/8 m 2 m\ - 5/4 m 2 s - 5/8 m 2 K s -l/8m A K + 3/4 s 2 ) 
+B(m 2 K ,m 2 K ,s) (1/4 m 2 s + 1/4 s - 1/4 s 2 ) 

+B(m 2 v , m 2 K , s) (-5/24 m 2 n m 2 K + l/8m 2 K s + 1/24 m A K ) + 5(m 2 , m 2 , s) (1/18 m 2 m^) 
+^(m 2 , m 2 K , s) (-3/4 m 2 + 1/2 m A n + 1/4 m A K ) 

+Bl(m% m 2 K , s) (1/4 m 2 m 2 K - 1/4 m 4 K ) . (B.4) 



+N; (s 2 ) + U 2 (8 S 2 ) + (1/(16tt 2 )) (1/6 s 2 ) 

+B(m 2 n , m 2 n , s) (1/2 m 2 - 3/4 m 2 s + 1/2 m 4 n - 1/4 m\ s + 1/4 s 2 ) 

+5(m 2 , m 2 K , s) (1/2 m 2 m 2 K - 1/2 m 2 s + 1/4 m% - 1/2 s + 1/4 + 1/4 s 2 ). 

(B.5) 
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Ml(s) = 

+B(ml, ml s ) (1/3 ml- 1/12 s) + B(ml, m 2 K ,s)(l/8 ml + 1/24 m\ - 1/24 s) 
+B(m 2 v , m 2 K , s) (-1/8 ml + 5/8 m\ - 1/8 s) + B[{ml, m 2 K , s) (-1/12 m 2 + 1/12 m 2 K ) 
+Bl{m% m 2 K , s) (1/12 m 2 - 1/12 m 2 K ). (B.6) 

Ml{s) = 0. (B.7) 

M 5 8 ( S )_= 

+B(ml, m 2 K , s) (-1/8 m 2 - 1/24 m\ + 1/24 s) + fi(m|, m|, s) (-1/3 m\ + 1/12 s) 
+B(m 2 , m 2 K , s) (1/8 m 2 - 5/8 m\ + 1/8 s) + ^(m 2 , m^, s) (1/12 m 2 - 1/12 m^) 
+Bl{m% m\, s) (-1/12 m 2 + 1/12 m 2 K ); (B.8) 

<(,) = -M|( S ). (B.9) 



K00 = 

+(N{ + iV 2 r ) (-2m>^ + 3m^-5mJ + m^- m^) 

(2 m 2 - 3mjs-m*- m\ s - m A K + 2 s 2 ) 
+iVJ(-8m>| f -2ra^ + 2m 4 I + 2m^-2m[) 
+7Vf(-4ra 2 m|- - 4 m 2 K s + 4 m 4 K ) 
+N^(-10mlm 2 K -2mls + 2ml + 4:m 2 K s-4m 4 K ) 
+iV 9 r (-2m^ + 2ra*+2ra^-2m 4 f( ) + (7V[ + 2N r n + N[ 2 ) (8 m 2 m 2 K ) 
+L[ (32 m 2 - 48 m 2 s + 80 m 4 n - 16 m 2 K s + 16 m^) 
+L2 (-16 ml m 2 K + 24 m 2 s + 8 + 8 m 2 K s + 8 m A K - 16 s 2 ) 
+L r 3 (16 m 2 - 24 m 2 s + 40 m% - 8 s + 8 m A K ) 
+LI (16 ml m 2 K + lQmls- 16 m*) + L\ (16 m 2 s - 16 m*) 
+ (1/(16tt 2 )) (11/9 m 2 ™^ + 4/9 m 2 s - 1/3 - 13/9 777^5 + 4/9 777^- 1/3 s 2 ) 
+A(ml) (-25/8 ml - 31/24 + 21/8 s) + A(m^) (-5/2 m 2 - 7/3 m| + 3s) 
+A(m 2 ) (-3/8 ml - 3/8 + 3/8 s) 

+B(ml, ml, s) (-ml m 2 K + 3/2 m 2 s - m\ + 1/2 m 2 K s - 1/2 s 2 ) 

+fi(m 2 , , s) (-m 2 m 2 K + m 2 s - 1/2 + s - 1/2 m A K - 1/2 s 2 ). (B. 10) 

K(«) = 

+ (iV[ + N r 2 + iV 3 r - 16L[ - 8L r 2 - 8L r 3 ) (s 2 ) 

+ (1/(16tt 2 )) (-l/18m 2 s + l/18m 2 K s + 1/6 s 2 ) 

+B(ml , ml, s) (5/4 m 2 s - 1/2 - 1/4 s - 3/4 s 2 ) 
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+B(ml, m 2 K , s) (-9/8 m 2 m\ + 3/4 m 2 s + 1/4 m\ + 1/8 m\ s + 3/8 m% - 1/2 s 2 ) 
+B(m 2 K , m 2 K , s) (-1/4 m 2 s - 1/4 s + 1/4 s 2 ) 

+5(m 2 ,m 2 oS ) (5/24 m 2 ™ 2 , - 1/8 m^s - 1/24 m^) + B(m 2 , m 2 , s) (- 1/18 m 2 ™ 2 ,) 
+Th(ml,m 2 K , s) (3/4 m 2 ™^ - 1/2 m 4 - 1/4 m^) 

+^(m 2 ,m 2 s) (-l/4m 2 m 2 , + 1/4 m^). (B.ll) 

M* 2 (s) = M!(s) = -M 9 8 ( S ). (B.12) 
The 27 plet ones depend as well on the quantity 

Pk* = -t^ J- (B.13) 



m" K — m 



M 27 (s) = 

+D[(-A/3mlm 2 K ) 

+D' 2 ((2/3 m 2 m\ + 10 m 2 s - 40/3 m*) + P/r^lO m\ s - 40/3 m°)) 
+D r 4 ((-11/2 m 2 + m 2 s + 7/3 + 3/2 m 2 K s - 5/6 m 4 K ) 

+ P K7T (-l0m 4 n s + 40/3mD) 
+D r 5 (29/6 ml m 2 K - 11 mis + 11 m 4 w - 3/2 s + 5/6 m^) 
+D r 6 (5/3 m 2 + 1/3 m 2 s - 1/3 - s) 

((-25/3 ml m 2 K - 3 m 2 s + 8 m£ + 6 m\ s - 10/3 m*.) 

+ PK W D?(-15m^ + 20m«)) 
+P^ 6 (-2/3 ml m 2 K + ml s - 2/3 m 4 + 1/3 s - 1/3 s 2 ) 
+Z^ 7 (4/3 m 2 - 2 m 2 s + 10/3 - 2/3 m\ s + 2/3 m 4 K - 4/3 s 2 ) 
+D r 28 (8 m 2 - 12 m 2 s + 12 m\ - 4 s + 4/3 m 4 K ) 
+D r 29 (-3/2 ml m 2 K + 9/4 m 2 s - 13/4 + 3/4 m 2 K s-7/12m 4 K + s 2 ) 
+D r 30 (3/2 ml m\ - 9/4 m 2 s + 1/4 - 3/4 m\ s - 5/12 m 4 K + 2 s 2 ) 
(4 m 2 m|V - 6 m 2 s + 4 m 4 - 2 s + 2 s 2 ) 
(-32 ml m 2 K + A8mls- 32 + 16 m 2 K s - 16 s 2 ) 
+L r 2 (16 m 2 - 24 m 2 s + 40 m 4 - 8 s + 8 - 16 s 2 ) 
+P 3 ((-58/3 m 2 m 2 K + 34 m 2 s - 1 18/3 m\ + 8 m 2 K s - 8 s 2 ) 

+ P Kt (40 m^s- 160/3 m^)) 
+U 4 (16 ml m 2 K -lQmls + 16 m*) 

+L£ ((8 ml m\ - 16 m 2 s + 68/3 m*) + P Kn (-20 s + 80/3 ml)) 

+ (1/(16tt 2 )) ((503/324 m 2 m 2 K + 197/216 m 2 s - 287/216 - 401/216 m 2 K s 

+ 401/648 m 4 K ) + P x?r (1/(16tt 2 )) (5/4 m 4 . s - 5/3 m°)) 
+A(ml ) ((-409/72 m 2 - 49/9 m 2 K + 29/3 s) + P Kn (-215/24 m 2 s + 215/18 m 4 ) 
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+ P 2 W (5/2mJs-10/3mJ)) 
+A{m 2 K ) ((-191/72 ml - 19/8 m\ + 37/6 s) + P^ (-325/24 m 2 s + 325/18 m%) 

+ ^(-5/2mJa + 10/3mJ)) 
+A(m 2 ) ((-13/12 m 2 - 5/8 + 3/2 s) + P Ku (-5/4 m 2 s + 5/3 m*)) 
+B(ml, ml, s) ((1/3 m 2 m\ - 13/6 m 2 s + 2/3 m\ - 5/6 m\ s + 3/2 s 2 ) 

+ P^ (-5/2 m 2 s 2 + 35/6 s - 10/3 m®)) 
+P(m 2 , m^, s) ((91/24 m 2 - 35/6 m 2 s + 10/3 - 25/12 m\ s - 7/24 

+ 19/8 s 2 ) + P Kn (15/8 m^ 2 -5mjs+ 10/3 m°)) 
+B(m 2 K , m 2 K , s) ((-13/12 m 2 s-7/8m 2 K s + 7/8 s 2 ) + P^ (5/8 m 2 s 2 - 5/6 m\ s)) 
+B{m 2 v ,m 2 K ,s) (5/24 m 2 m^ - 1/8 m\ s - 1/24 m%) + B{m 2 , m 2 , s) (-l/18m 2 m^) 
+P^(m 2 , m^, s) (-4/3 m 2 + 1/3 - 2/3 m%) 

^B\{m\,m 2 K , s) (-1/4 mlm 2 K + 1/4 m^). (B.14) 
Ml\s) = 

+D r 27 (-2/3 s 2 ) + P 2 r 8 (-4/3 s 2 ) + D r 29 (7/12 s 2 ) + D r 30 (5/12 s 2 ) + (-8 s 2 ) 
+P K7T L r 3 (10/3 m 2 s 2 ) 

+ (1/(16tt 2 )) ((-55/648 m 2 s + 55/648 s + 8/27 s 2 ) + P K7r (-5/24 m 2 s 2 )) 
+P(m 2 , ml, s) ((-1/12 m 2 - 11/12 m 2 s + 1/3 m% + 1/24 s + 3/8 s 2 ) 

+ P Kw B(ml, ml, s) (-5/8 m 2 s 2 + 25/12 s - 5/3 m£)) 
+B{m 2 n ,m 2 K ,s) ((13/4 m 2 ™^ -67/16 m 2 s + 11/6 m\ -217/96m 2 K s + 1/6 

+ 67/32 s 2 ) + P K7T B(ml, m\, s) (25/32 m 2 s 2 - 5/3 m\ s + 5/6 m®)) 
+P(mJ, m^, s) ((5/72 m 2 - 5/12 m 2 s + 5/18 - 5/32 s - 5/72 

+ 5/32 s 2 ) + P K7T (5/32 ml s 2 - 5/12 m\ s + 5/18 m£)) 
+B[{ml,m 2 K ,s) (-15/16 m 2 + 5/24 - 25/48 m^) 

+P^(m 2 ,,774,s) (5/16m 2 m^ + 5/48m^). (B.15) 
Mf{s)_= 

+B(ml, ml, s) ((1/2 m 2 - 1/8 s) + P K * (5/24 m 2 s - 5/6 ml)) 

+B(m 2 n , m\, s) (67/144 m 2 + 43/288 - 43/288 s + P^ (-5/288 m 2 s + 5/72 m*)) 
+P(m^, m\, s) ((-5/18 m 2 + 5/9 m\ - 5/36 s) + P KlT (5/72 m 2 s - 5/18 m*)) 
+P(m 2 , m^, s) ((-7/24 m 2 + 5/32 - 1/32 s) + P Kn (5/96 m 2 s - 5/24 m*)) 
+Bl(ml,m 2 K ,s) (-19/72 m 2 +43/1447/1^) +B~ 1 (m 2 r) ,m 2 K ,s) (1/18 m 2 - 1/48 m\). 

(B.16) 

Mf ( a ) = 
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+D r 2 (10/3 mls) + P K7T (10/3 m*s)) + D r 4 (-10 ra 2 s- 5/2 s) + P K * (-10/3 m*s)) 
+P£ (20/3 m 2 s + 5/2 s) + P>£ ((-25/3 m 2 s - 10 s) + P KlT (-5 s)) 
+D r 28 (10 m 2 s + 10/3 m^s- 10/3 s 2 ) + (P^ - Py (-5/4 m 2 s - 5/12 77^5 + 5/12 s 2 ) 
+P, ((10/3 m 2 s) + P x?r (-10/3 m 2 s 2 + 40/3 s)) 
+L\ (64/9 m 2 m 2 ,) + L r 5 (P K7T (-20/3 s)) 

+ (1/(16tt 2 )) (5/18m 2 s + 205/108 m 2 K s + 55/54 s 2 ) + P x?r (5/24m 2 s 2 + 5/12m^s)) 
+A(m 2 ) ((-1105/72 s) + P Kw (-95/24 m 2 s) + P 2 Kn (5/6 s)) 
+A(m 2 Y ) ((-335/36 s) + Px. (-85/24 m 2 s) + P 2 , (-5/6 m\ s)) 
+A(m 2 ) ((-35/24 s) + P Kw (-5/12 m 2 s)) 

+P(m 2 , m 2 , s) ((5/4 m 2 - 35/6 m 2 s + 25/6 mJ-5/8m^ + 15/8 s 2 ) 

+ P^ (5/8 m 2 s 2 - 25/12 s + 5/3 ml)) 
+B(ml, m 2 K , s) ((-5/12 m 2 - 5/16 m 2 s + 5/12 - 55/96 m 2 K s + 5/12 

+ 5/32 s 2 ) + P Kw (-25/32 m 2 s 2 + 5/3 ml s - 5/6 ml)) 
+B(m 2 , m 2 K , s) (-5/72 m 2 + 5/12 m 2 s - 5/18 + 5/32 m 2 K s + 5/72 m A K 

- 5/32 s 2 ) + P Kn (-5/32 m 2 s 2 + 5/12 m\ s - 5/18 ml)) 
+B^(ml,m 2 K , s) (15/16 m 2 - 5/24 m* + 25/48 m^) 

+P^(m 2 ,m^,s) (-5/16 m 2 - 5/48 m A K ). (B.17) 
Mf( S ) = 

+P 2 r ((10/3 ml) + P X7r (10/3 ml)) + D\ ((-10 m 2 - 5/2 m\) + P^ (-10/3 ml)) 

+D r 5 (20/3 m 2 + 5/2 m\) + ((-25/3 m 2 - 10 m 2 ,) + P KlT (-5 m*)) 

+P 2 r 8 (10/3 s) + P 2 r 9 (-5/12 s) + D r 30 (5/12 s) + P K „ (L r 3 - 2I/ 5 ) (10/3 m 2 s) 

+ (1/(16tt 2 )) ((280/81 ml + 1945/648 m 2 K - 55/54 s) + P x?r (-5/24 m 2 s + 5/4 m*)) 

+Z(m 2 ) ((-1105/72) + P Kw (-95/24 m 2 ) + P^ (5/6 m*)) 

+Z(m 2 ,) ((-335/36) + P^ (-85/24 m 2 ) + P 2 w (-5/6 m*)) 

+Z(m 2 ) ((-35/24) + P Kn (-5/12 m 2 )) 

+P(m 2 , m 2 , s) (-5/6 m 2 + 5/24 s) + P Kw (5/72 m 2 s - 5/18 m*)) 
+P(m 2 , m 2 K , s) (-67/144 m 2 - 43/288 m 2 K + 43/288 s + P Kn (5/288 m 2 s - 5/72 ml)) 
+B(m 2 K , m 2 K , s) ((-5/18 m 2 - 7/9 + 7/36 s) + P Ku (5/72 m 2 s - 5/18 m*)) 
+B(m 2 , m 2 K , s) ((7/24 m 2 - 5/32 m\ + 1/32 s) + P K7T (-5/96 m 2 s + 5/24 m*)) 
+PT(m 2 , m 2 K , s) (19/72 m 2 - 43/144 m\) + PT(mJ, m^, s) (-1/18 m 2 + 1/48 m\). 

(B.18) 

M 27 ( S |= 

+P(m 2 , ml , s) ((-19/18 m 2 + 19/72 s) + P x?r (5/24 m 2 s - 5/6 ml)) 
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+B(m 2 n , m 2 K , s) (-67/144 m 2 - 43/288 m\ + 43/288 s + P X7r (5/288 m 2 s - 5/72 m*)) 
+B(m 2 K , m 2 K , s) ((-5/9 m 2 - 5/9 m\ + 5/36 s) + P Ku (5/36 m 2 s - 5/9 m*)) 
+P(tt7 2 , m^, s) ((7/24 m 2 - 5/32 + 1/32 s) + P KlT (-5/96 m 2 s + 5/24 m*)) 
+P7~(777 2 , m^, s) (19/72 m 2 - 43/144 m^) + PTYm 2 , m^, s) (-1/18 ml + 1/48 m^). 

(B.19) 

Mll{s) = 

(D[ + 2D r 2 ) (8/3 m 2 m^) + D r 4 (U/3 mlm\ - 47/3 m 2 s + 47/3 - 13/3 m\ s + m 4 K ) 

+D r 5 (-10 ml m\ + 47/3 m\ s - 47/3 m\ + 13/3 s - m A K ) 

+D r 6 (-2/3 m 2 m\ + 1/3 m 2 s - 1/3 m\ - m\ s + m A K ) 

+D r 7 (50/3 ml m\ - 34/3 mls + 34/3 m\ - 52/3 s + 4 m^) 

+P£ 6 (-2/3 mlm 2 K + mls- 5/3 + l/3ra^-l /3 m^) 

+£>5 7 (4/3 m>^-2m^- 2/3 m 4 n - 2/3 s - 2/3 + 4/3 s 2 ) 

+D r 28 (-12 m 2 m\ + 18 m 2 s - 18 m\ + 6 m\ s - 2 - 4 s 2 ) 

+P£ 9 (8/3 m 2 m*. - 4 m 2 s + 8/3 - 4/3 m\ s + 4/3 s 2 ) 

+D r 30 (-8/3 m 2 + 4 m 2 s - 20/3 + 4/3 m\ s - 4/3 m A K ) 

+D r 31 (-8/3 m 2 m\ + 4 m 2 s - 20/3 m\ + 4/3 m 2 K s- 4/3 m^) 

(64/3 ml m 2 K -32mls + 160/3 m£ - 32/3 m\ s + 32/3 m%) 
+L r 2 (-32/3 m 2 + 16 m 2 s + 16/3 m\_ + 16/3 s + 16/3 - 32/3 s 2 ) 
+L r 3 (32/3 m 2 m\ - 16 m 2 s + 80/3 - 16/3 m 2 K s + 16/3 m^) 
+LJ (32/3 ml m\ + 32/3 m 2 s - 32/3 m£) + L\ (32/3 m 2 s - 32/3 m\) 
+ (1/(16tt 2 )) (-643/162 m 2 - 49/324m 2 s + 31/108m^ + 1273/324777^5 

-439/324 777^ + 34/27 s 2 ) 
+A(ml ) ((2075/72 m 2 + 493/72 m 2 K - 329/12 s) + P Kn (-35/24 m 2 s + 25/9 m*)) 
+^(777^) ((575/36 m 2 + 149/36 - 49/3 s) + P Ku (5/6 m 2 s - 10/9 m*)) 
+A(m 2 v ) ((29/24 m 2 + 19/24™^ - 9/4 s) + (5/8 m 2 s - 5/3mJ)) 
+B(ml, ml, s) {ml m\ - 13/2 m 2 s + 13/3 m£ - 1/2 s + 13/6 s 2 ) 
+B(ml, m 2 K , s) (8/3 m 2 - 13/3 m 2 s + 13/6 - 8/3 m\ s + 1/2 m 4 K + 13/6 s 2 ). 

(B.20) 

M 2 l(s) = 

+ (D r 26 + 2D r 27 + QD 28 + AD r 30 + 4D r 31 ) (1/3 S 2 ) + (2L[ + U 2 + L\) (-16/3 S 2 ) 
+ (1/(16tt 2 )) (-11/162 777 2 s + 11/162 777^5- 17/27 s 2 ) 
+B(m 2 n , ml, s) ((-5/12 m 2 m\ + 5/2 m 2 s — 2 ra„. — 1/6 777^5 - 1/2 S 2 ) 
+ P^ (-5/4 777 2 s 2 + 35/12 777^ s - 5/3 ml)) 
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+B(ml , m 2 K , s) ((7/24 m 2 m\ - 13/24 m 2 s + 7/12 m\ + 19/48 m\ s - 3/8 m% 

- 1 /48 s 2 ) + P K7T (15/16 m 2 s 2 - 5/2 m\ s + 5/3 ml)) 
+B(m 2 K , m 2 K , s) (-3/8 m 2 s - 3/8 s + 3/8 s 2 ) 

+B(m 2 v m 2 K , s) ((-5/72 m^m^ - 5/6 m 2 s + 5/9 - 3/16 m^s - 7/72 

+ 5/16 s 2 ) + P x?r (5/16 ml s 2 - 5/6 s + 5/9 m£)) 
+P(m 2 , m 2 , s) ((-13/36 m 2 + 5/12 m 2 s - 5/9 m*) + P^ (5/12 s - 5/9 m£)) 
+P^(m 2 , s) (-13/24 m 2 m 2 K + 1/12 - 3/8 m 4 K ) 

+B\{m%m 2 Kl s) (7/8mlm 2 K - 1/24 m%). (B.21) 



M 2 J( S |= 

+P(m 2 , ml , s) (-13/9 m 2 + 13/36 s) 

+P(m 2 , m^, s) ((-8/9 m 2 - 41/144 m\ + 41/144 s) + P Kn (5/144 m 2 s - 5/36 m*)) 
+B(m 2 K , m^, s) ((5/9 m 2 - 5/18 rr& + 5/72 s) + P K7r (-5/36 m 2 s + 5/9 m*)) 
+P(m 2 / , m^, s) ((-7/24 m 2 - 55/48 m|- + 11/48 s) + P K7T (5/48 m 2 s - 5/12 m*)) 
+P^(m 2 , m^, s) (1/2 ml - 41/72 m|) + Bl{m% m 2 K , s) (-1/12 m 2 + 11/72 m\). 

(B.22) 
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